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^ ■ We perform a comprehensive and detailed comparison of the physics reach of Beta-beam neutrino 
! experiments between two pairs of plausible source ions, ( 8 B, 8 Li) and ( 18 Ne, 6 He). We study the 
optimal choices for the baseline, boost factor, and luminosity. We take a 50 kton iron calorimeter, 
a la ICAL@INO, as the far detector. We follow two complementary approaches for our study: 
(i) Fixing the number of useful ion decays and boost factor of the beam, and optimizing for the 
sensitivity reach between the two pairs of ions as a function of the baseline, (ii) Matching the 
shape of the spectrum between the two pairs of ions, and studying the requirements for baseline, 
boost factor, and luminosity. We find that for each pair of ions there are two baselines with 
very good sensitivity reaches: a short baseline with L [km]/7 ~ 2.6 ( 8 B+ 8 Li) and L [km]/7 ~ 0.8 
( 18 Ne+ 6 He), and a long "magic" baseline. For 7 ~ 500, one would optimally use 18 Ne and 6 He 
at the short baseline for CP violation, 8 B and 8 Li at the magic baseline for the mass hierarchy, 
and either 18 Ne and 6 He at the short baseline or 8 B and 8 Li at the magic baseline for the sin 2 2# 13 
discovery. 
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1 Introduction 



Spectacular results from a series of neutrino oscillation experiments over the last four decades [1, 
2, 3, 4, 5, 6, 7] have paved the way for the "golden" age of neutrino physics. Determining the 
hitherto unknown mixing angle #13, the CP phase Sep, and the sign of Am^, i.e., the neutrino 
mass hierarchy 1 , have emerged as the next frontier in this field. All these three quantities can be 
probed by experimentally measuring the so-called "golden" channel [8] oscillation probability P efJi 
(or its T-conjugate channel P^ e )- A series of ambitious projects are under discussion which plan 
to use this oscillation channel. The on-going and near future experiments include the MINOS 
experiment in the US [6], and the CNGS experiments ICARUS [9] and OPERA [10] in Europe. 
Next experiments in line will be T2K in Japan [11] and NOz/A in US [12]. All these experiments 
will use muon neutrino beams from conventional accelerator sources in order to observe P^, e . 
Collectively and in combination with short-baseline reactor experiments, such as Double Chooz 
[13], these experiments are expected to improve the bound on #13 to about sin 2 2^13 < 0.01 (90% 
CL) [14]. The mass hierarchy and CP violation, though in principle accessible using the combined 
data from the T2K and NOz/A experiments, can be determined only for values of sin 2 2#i 3 (true) 
close to the current bound and for some fraction of the possible values of the CP phase <5cp (true). 2 
The sensitivity of these experiments is mainly restricted by statistics, while for larger luminosity 
set-ups, the intrinsic v e background poses a natural limitation for experiments sensitive to z/ M 
oscillations into v e . Therefore, if Nature has not been very kind we will need larger experiments 
to complete our understanding of the neutrinos, possibly using an alternate technology. 

In order to access small values of sin 2 2# 13 , there are several requirements: One needs to 
have low backgrounds, more statistics, and reduced systematical uncertainties. As far as the low 
background requirement is concerned, it is an advantage to use a pure flavor neutrino beam without 
any intrinsic beam contamination. One such approach has been proposed by Piero Zucchelli [15]: 
Radioactive nuclides are created by impinging a target by accelerated protons. These unstable 
nuclides are collected, fully ionized, bunched, accelerated and then stored in a decay ring (see 
for e.g. [16, 17]). The decay of these highly boosted ions in the straight sections of the decay 
ring produces the so-called Beta-beam. An alternative approach is the so-called Neutrino Factory 
(NuFact) [18]. It involves producing, collecting, cooling, accelerating, and circulating muon packets 
in a storage ring. The decay of accelerated muons (antimuons) in the straight sections of the storage 
rings produce and v e [p^ and z/ e ) beams. The presence of the u e (z/ e ) in the beam allows for the 
observation of the Pep, {Pen) oscillation probability in the far detector. Since both z/ M (or z/ M ) from 
the original beam as well as z/ M (or z/^) from the oscillated v e (or z/ e ) will be arriving at the detector, 
it must have the ability to distinguish one from the other. The most accepted candidate is the 
magnetized iron detector, though there are several proposals with more expensive and elaborate 
designs and therefore better performance [19]. Statistics can be increased by a higher beam power 
and the size and efficiency of the detector. Beam-related systematic uncertainties can be reduced 

1 Though we call this the neutrino mass hierarchy, what we mean is basically the ordering of the neutrino mass 
states. Therefore, our discussions are valid for both hierarchical as well as quasi-degenerate mass spectra. We 
define Am| = mf —rrij and refer to s<7n(Am§ 1 ) as the neutrino mass hierarchy - s<7n(Am| 1 ) > is called "normal 
hierarchy" (NH) while sgn^m^) < is called "inverted hierarchy" (IH) . 

2 We distinguish between the "true" values of the oscillation parameters, which are the values chosen by Nature, 
and their fitted values. Throughout this paper we denote the true value of a parameter by putting "(true)" after 
the symbol for the parameter. 
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to a large extent by working with a two detector set-up, one very close to the beam line and another 
serving as the far detector. The systematic uncertainties coming from the lack of knowledge of 
the neutrino-nucleus interaction cross-sections are another important source of error. These can 
be controlled to some degree by the near-far two detector set-up, but they cannot be canceled 
completely [20]. Beam- related backgrounds are extremely small for the NuFact and Beta-beam 
experiments because they either use leptonic decays (NuFact) or a flavor-pure beam (Beta-beam). 
The detector backgrounds coming mainly from neutral current interactions and mis-identification 
of particles, can be reduced by imposing intelligent cuts. The atmospheric neutrino backgrounds, 
which can be important for Beta-beams at lower energies, can be suppressed using timing and 
directional information. 

The performance and physics reach of these expensive and ambitious experiments have been 
the subject of much discussion for the past few years [21]. There has been a plethora of papers on 
this issue 3 , most of which have addressed the problem of "parameter degeneracies". Even if both 
neutrinos and antineutrinos are used, there are three types of discrete degeneracies in the golden 
channel: 

1. the (#i3,<5 C p) intrinsic degeneracy [23], 

2. the (sgn(Am| 1 ), 5cp) degeneracy [24], 

3. the (#23, 7r/2 — #23) degeneracy [25]. 

Together they can result in up to eight-fold degenerate solutions [26], severely deteriorating the 
sensitivity of the experiment. The variety of suggestions to solve this problem includes combining 
data from several experiments observing the golden channel, but with different baselines L and 
neutrino energies E [23, 27, 28], combining data from accelerator experiments observing different 
oscillation channels [29, 30, 31], combining the golden channel data with those from atmospheric 
neutrino [32, 33] or reactor antineutrino experiments [34]. A particularly attractive way of com- 
pletely resolving at least two of the three degeneracies is to perform the experiment at the "magic 
baseline" [35, 36, 37]. This magic baseline reflects the characteristic oscillation wavelength corre- 
sponding to Earth matter. One can show that at this baseline, for reasonably small values of #13, 
the ^cp-dependent terms vanish. The <5cp-dependence is therefore reduced and one can eliminate 
the (6*i3,5cp) and (sg^AnigJ, 5 CP ) degeneracies, resulting in tremendous sensitivity to 9 13 and 
sgn(Amg 1 ). The sensitivity reach of a NuFact experiment at the magic baseline can be found in 
[35, 38]. The idea for a magic baseline Beta-beam experiment with a similar performance for #13 
and sgn(Am| 1 ) was put forth in Refs. [39, 40, 41]. 

The sensitivity reach of an experiment depends crucially on beam, baseline, and detector 
properties. It is therefore important to ask which beam, baseline and detector set-up would 
qualify as the optimal choice in order to obtain the best sin 2 2#i3 reaches for the three quantities 
we have set out to measure, i.e., #13, sgn(Am| 1 ), and <5cp. For the NuFact, this detailed exercise 
was performed in Ref. [38]. The sensitivity to each of the three parameters mentioned above was 
studied as a function of the baseline and muon neutrino energy. It was demonstrated that the 

3 A summary of the potential of selected NuFact and Beta-beam set-ups have been compiled by the physics 
working group of the International Scoping Study for a future Neutrino Factory, Superbeam and Beta-beam, in 
their report [22]. 
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minimal muon neutrino energy acceptable for the magnetized iron detector was about 20 GeV. 
The optimal baseline for probing 13 and sgn(Am| 1 ) is the magic baseline, whereas best sensitivity 
to CP violation is expected at L ~ 3000 to 5000 km [38]. 

For Beta-beams, a variety of plausible set-ups have been proposed in the literature [33, 39, 
40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53]. The proposal which poses minimal challenge 
for the Beta-beam design, is the commonly called CERN-MEMPHYS project [33, 42, 43]. It 
proposes to use the EURISOL ion source to produce the radioactive source ions 18 Ne and 6 He, 
and demands a Lorentz boost factor 7 ~ 100 for them, which can be produced using the existing 
accelerator facilities at CERN. The far detector MEMPHYS, a megaton water detector with 
fiducial mass of 440 kton, will have to be built in the Frejus tunnel, at a distance of 130 km from 
CERN. Another possible Beta-beam set-up using water detector but higher boost factors and an 
intermediate baseline option was put forth in [46, 47] (see also [51]). In these papers authors 
have used a high 7 18 Ne and 6 He Beta-beam option at CERN and 440 kton fiducial volume 
water detector at GranSasso or Canfranc, which corresponds to L — 730 and 650 km respectively. 
Excellent sensitivity to #13 and CP violation is expected [46, 47] from this proposal. Another 
high performance set-up proposed in [39, 40, 41] would use a high 7 Beta-beam and a magnetized 
iron detector placed at a distance close to the magic baseline. Since high energy neutrinos are 
mandatory for achieving near-resonant matter effects required for the desired performance of this 
set-up, one needs to employ an alternative set of source ions, 8 B and 8 Li [54, 55]. The end-point 
energy of these ions are larger than those of 18 Ne and 6 He by factor of about 3.5, and hence optimal 
neutrino energies can be obtained with a 7 between 350 and 650. As one possible option, the Beta- 
beam could be targeted from CERN towards the India-Based Neutrino Observatory (INO) [56]. 
The CERN to INO distance corresponds to 7152 km, which is almost magic. Therefore, this 
experiment yields sensitivity to #13 and sg^Am^), both of which could be outperformed only by 
the NuFact experiment at the magic baseline distance. Set-ups with a neutrino beam from CERN 
to GranSasso or CanFranc [44] , from CERN to Boulby mine [53] and from Fermilab [48] (L ~ 300 
km) have also been proposed, and their sensitivity reach has been explored. Set-up with two sets 
of source ions with different boost factor for each set but with the same baseline was proposed 
in [51]. In [45] the authors consider the complementary situation where they take only one set of 
source ions, 8 B and 8 Li, with 7 = 350 and two different baselines, L = 2000 km and 7000 km. 
Very high gamma Beta-beam options have been studied in Refs. [39, 47, 49]. The physics potential 
of low energy Beta-beam option was probed in [57]. A comparison of the physics reach among 
different Beta-beam experimental proposals can be found in Ref. [50, 58]. 

In contrast to earlier works, we study the baseline optimization as a function of the ion pair 
used 18 Ne+ 6 He or 8 B+ 8 Li, and we discuss the impact of the luminosity. In addition, we perform 
a simultaneous optimization of L and 7. All comparisons are performed for the same detector, 
which is a 50 kton magnetized iron calorimeter. Note that the magnetization of the detector, 
which is mandatory for the Neutrino Factory, is only used for a reduction of the backgrounds. 
We study the performance with respect to 9i 3 , sg^Am^), and CP violation, paper is organized 
as follows. We describe the Beta-beam experiment and our analysis procedure in Section 2. In 
Section 3, we optimize the baseline L of the experiment for fixed sets of 7's. In Section 4, we 
then perform a simultaneous optimization over L and 7, and we discuss the requirements for a 
similar beam spectrum. In Section 5, we then show the impact of the Beta-beam luminosity, and 
determine the best combination of L, 7, and N@ (the number of useful ion decays per year). Our 
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Ion 


r(s) 


E (MeV) 


/ 


Decay fraction 


Beam 


18 Ne 


2.41 


3.92 


820.37 


92.1% 


v e 


tRe 


1.17 


4.02 


934.53 


100% 


Ve 


1* 


1.11 


14.43 


600872.07 


100% 


Ve 


|Li 


1.20 


13.47 


425355.16 


100% 





Table 1: Beta decay parameters: lifetime r, electron total end-point energy E , /-value and decay 
fraction for various ions [60]. 



conclusions can be found in Section 6. 

2 Simulation of the Beta-beam Experiment 

Here we describe the experimental set-up for our proposed Beta-beam facility. We give the details 
of the flux and detector set-up we have used in our analysis. 

2.1 The Flux 

A Beta-beam [15] is an intense and highly collimated source of pure v e or v e flux, produced from 
the decay of beta unstable radioactive ions. These unstable ions are created by impinging a target 
with high energy protons. Subsequently, the ions are collected, bunched, accelerated and stored 
in a decay ring. The standard design of the decay ring comprises of a racetrack shaped tunnel. 
When the ions decay along the straight sections, they produce a v e or v e beam. 

This neutrino beam would be very suitable for precision experiments because it is mono-flavor 
and hence beam related backgrounds are almost absent. The neutrino spectrum depends only on 
the beta decay total end-point energy E and the Lorentz boost of the radioactive ions 7. The 
spectral shape can therefore be very well determined. The flux normalization is given directly 
by N/3, the number of useful ion decays per year in the straight section of the storage ring. The 
standard numbers taken for the 18 Ne and 6 He are 1.1 x 10 18 (is e ) and 2.9 x 10 18 (P e ) useful decays 
per year, respectively [59]. Wherever not explicitly mentioned, these reference numbers of useful 
ion decays for u e and v e are chosen. Note, however, that new ideas suggest luminosities higher even 
by a factor of ten or so, depending on the isotopes used, by using a recirculating ring to improve 
the performance of the ion source [54, 55]. The total luminosity is given by the product of useful 
ion decays/year x running time x detector mass x detection efficiency Throughout the study 
we will consider five years of neutrino and five years of antineutrino running. Since the reference 
luminosity might not be reachable for different reasons, or it may be much higher because of a 
better ion source, larger detector, etc., we will include it as a parameter in this study. 

The beam divergence is controlled by the Lorentz boost 7. Hence by increasing 7, we can 
produce a higher beam collimation and increase the beam intensity along the forward direction 
oc 7 2 . However, note that though the intensity can be increased by choosing either a higher Np or 
7, they might produce very different neutrino beams. While increasing Np merely increases the 
overall normalization of the flux, increasing 7 increases both the total flux as well as the average 



5 




2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 



Energy (GeV) Energy (GeV) 




Energy (GeV) Energy (GeV) 



Figure 1: The unoscillated Beta-beam flux spectrum arriving at a detector placed at the magic 
baseline. The upper panels are for 8 B (left panel) and 8 Li (right panel), while the lower panels are 
for 18 Ne (left panel) and 6 He (right panel). 



energy of the beam. This dependence of the beam flux on 7 is illustrated in Fig. 1. Notice that 
we assume the same 7 for both neutrino as well as antineutrino modes. 

Another crucial aspect associated with Beta-beams is the choice of the beta unstable ion. The 
properties that a suitable ion should have include a high production yield, large decay fraction, 
reasonably long lifetime, and preferably lower Z/A ratio. The most widely discussed ions are 
18 Ne and 6 He, which would produce a v e and v e beam respectively. The details of the source 
ions for Beta-beams are given in Table 1. The other pair of ions which have been proposed as an 
alternative to 18 Ne and 6 He, are 8 B and 8 Li [54, 55]. The main difference between 8 B+ 8 Li compared 
to 18 Ne+ 6 He is the higher end-point energy (see Table 1). The factor ~ 3.68 (3.35) difference in 
end-point energy ensures that for the same peak (anti)neutrino energy, approximately given by 
7-E'o, the 7 required for 8 B ( 8 Li ) will be 3.68 (3.35) times smaller than that needed for 18 Ne ( 6 He). 
Since we assume the same 7 for both ions within each pair, we use the average difference in the 
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total end-point energy (3.35 + 3.68)/2 ~ 3.5 to estimate the effects of 7. The isotope dependence 
can be also seen in Fig. 1. The figure shows that the peak energy is approximately given by jE . 

Let us discuss now the conditions to obtain a similar spectrum (including normalization) when 
using different isotope pairs at the same baseline. The purpose of this exercise is to derive the 
conditions under which we produce matching neutrino energies and fluxes, and therefore deal with 
the same physics (including neutrino energies, statistics, and, especially, matter effects). If we 
neglect effects of the endpoint in the beta beam spectrum (i.e., E ^> m e ), we know from the 
beta beam flux formula that the peak energy of the spectrum is approximately given by jEq, and 
the total flux is proportional to Np^ 2 . In order to obtain a spectrum with the same peak energy 
and normalization for two different isotopes with very different endpoint energies (such as 8 B and 
18 Ne), we therefore have the following conditions (modulo endpoint effects): 

jvf ~ \E<V ) ' 7« 2 > " £f 

From these matching conditions, one obtains the ratio of /N^ and j^ 2 ^ needed for the 
source ions. Therefore, using an isotope with a higher endpoint energy allows for a lower 7 to 
obtain the same neutrino energies. However, in order to get the same flux, the useful isotope 
decays have to be adjusted quadratically. For our pairs of isotopes, we have E® +Ll ~ 3.5 • E$ e+He . 
Therefore, from Eq. (1) we obtain the conditions 

Np +Li ~ 12 • Np e+He , ^ Ne + He ~ 3.5 . 7 B + Li 

in order to obtain the same neutrino flux spectrum. Note that the number of useful ion decays Np 
represents, to first approximation, an ion source degree of freedom, whereas the 7 is an accelerator 
degree of freedom. 4 Each can be adjusted with completely different technical challenges. The actual 
optimization between higher 7 versus higher isotope production rates depends on individual cost 
and machine aspects, and cannot be done at this place [61]. 

The optimal baseline depends crucially on the choice of source ions and the boost factor. 
For shorter baselines one is away from the matter resonance, and hence the flux arriving at the 
detector is proportional to 1/L 2 . If one wants to stay at the oscillation maximum in vacuum, one 
has L/E = const., and therefore L oc 7. Since the cross sections are proportional to ~ E oc 7 
for deep inelastic scattering (DIS) processes, one has an overall 1/L 2 X7X7 2 = 7 scaling of the 
event rates in the DIS regime. Close to matter resonance, the flux at the detector hardly falls as 
a function of L, which means that longer baselines might be preferred. This qualitative discussion 
of the baseline dependence does not take into account the non-trivial dependence of the oscillation 
probabilities on the oscillation parameters, and the intrinsic degeneracies. It is the purpose of this 
work to study this dependence. 

From the discussion above we see that one has to optimize for the Beta-beam flux itself by a 
judicious choice of: 

• The types of ions and their end-point energy. 

4 Thcrc is, however, a non-negligible effect of Np on the accelerator by increasing the number of ions per bunch 
(or the number of bunches). 
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• Lorentz boost factor, 7. 

• Number of useful ion decays per year, Np. 

It is clear that every choice of ion, 7, Np and L will give a different physics reach for the experiment. 
The choice of 7 and Np will determine the initial Beta-beam flux for a given choice of the source 
ions. This is what we would call the "input" of the experiment. What finally determines the 
physics reach of the experiment is the number of events seen in the detector, and the potential to 
resolve correlations and degeneracies. We will call this the "output" of the experiment. The aim 
of course is to maximize the "output" . However, there are practical limitations on stretching the 
"input" possibilities. Keeping these in mind, we study the comparative sensitivity reach of the 
greenfield Beta-beam set-ups in two ways: 

1. By fixing the input and comparing the output. 

2. By fixing the output and comparing the required input. 

We will use approach (1) in section 3, where we make a comparison between the sensitivity reach 
of the experiment using either the 8 B and 8 Li combination or the 18 Ne and 6 He combination, as a 
function of L. Sections 4 and 5 are more in the spirit of approach (2). 

2.2 The Detector 

We are interested in measuring the golden channel probability P efl . Since we have a v e (p e ) flux in 
the beam, we need a detector which is sensitive to muons (antimuons). The detector should have 
a suitable energy threshold, depending on the energy spectrum of the Beta-beam. In addition, 
it should have a good energy resolution and low backgrounds. There are a number of detector 
technologies that have been considered in the literature. For the low energy Beta-beams, water 
Cerenkov detectors are the most widely chosen, mainly because of their low energy threshold and 
large size. This is a very well known and tested detector technology. In addition, the detector 
can be relatively easily upgraded; typically megaton-sizes [62, 63, 64] appear in the literature. 
However, the backgrounds in this detector are generally larger than in other detector types. For an 
intermediate 7 Beta-beam, a Totally Active Scintillator Detector (TASD) is a possible technology. 
This is the option chosen and studied by the NOz/A collaboration [12]. The third kind of detector 
technology, which has been studied extensively and which is currently being used by MINOS, 
is the magnetized iron calorimeter. A larger version is envisaged to come up soon at the INO 
facility in India [56]. In this paper it will be referred to as ICAL@INO. For both TASD and 
magnetized iron detectors the background rejection is typically considered to be better than for 
water Cerenkov detectors. 

In this paper we use, for the sake of simplicity, only one type of detector for both types of ions 
and all values of 7. We use an ICAL@INO type of detector configuration [56]. We give the details 
of our detector specifications in Table 2. The charge identification efficiency is incorporated since 
that helps in reducing the neutral current backgrounds. The number of (anti)muon events in the 
detector is given by 

„ 77 ,. 77 max 

pmax f A 

N l = Tn n f ID e dE * dE A cj>(E) a^E) R(E,E A ) P eix {E) , (3) 

JO JE™™ 
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1 PpV 

1 VorC V 


ueteccion H/inciciicy ^ej 


oU/o 


Charge Identification Efficiency (f/n) 


95% 


Detector Energy Resolution function (er) 


0.15E 


Bin Size 


I GeV 


NC Background Rejection 


0.0001 


Signal error 


2.5% 


Background error 


5% 



Table 2: Detector characteristics for neutrinos/antinuetrinos used in the simulations. The bin size 
is kept fixed, while the number of bins is varied according to the maximum energy. 



where T is the total running time (taken as five years), <f)(E) is the unoscillated flux at the detector, 
e is the detector efficiency, n n are the number of target nucleons in the detector, fj^ is the charge 
identification efficiency and R(E, Ea) is the energy resolution function of the detector, for which 
we assume a Gaussian function. For muon (antimuon) events, a v is the neutrino (antineutrino) 
interaction cross-section. The quantities E and Ea are the true and reconstructed (anti) neutrino 
energy respectively. 



2.3 Neutrino Propagation and Simulation Details 

The expression for P efJi in matter [65, 66, 67], up to second order terms in the small quantities #13 
and a = Am^/Am^, is given by [8, 68] 

,„ 2 „ sin 2 [(l-i)A] 
P eil ~ sin 2 ^3 sin 2 2^3- 



where 



± a sin 2#i3 sin 2#i 2 sin 2# 23 sin <5cp sin( A) 



+ a sin 29 13 sin 29 12 sin 2# 23 cos 5qp cos(A) 



{l-Af 

sin(iA) sin[(l - A) A] 



A (1-A) 
sin(iA) sin[(l - A) A] 



A (1-A) 

+ a 2 cos 2 g 23 sin 2 2g 12 Sin2( ^ A) , (4) 

A 2 



a - Am l L A=— (5) 
A= AE ' A ~ Am^ (5) 
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lAm^true)! = 2.5 x 10" 3 eV 2 


a(Am 2 31 ) = 


1.5% 


sin 2 2^ 2 3 (true) | = 1.0 


<r(sin 2 29 23 ) 


= 1% 


Am 2 x (true) = 8.0 x 10" 5 eV 2 


a(Am 2 l ) = 


2% 


sin 2 # 12 (true) = 0.31 


a(sin 2 12 ) = 


= 6% 


p(true) = 1 (PREM) 


a(p) = 5% 





Table 3: Chosen benchmark values of oscillation parameters and their la estimated errors. The 
last row gives the corresponding values for the Earth matter density. 

and A = ±2\/2G F N e E is the matter potential, given in terms of the electron density N e and 
(anti)neutrino energy E; the plus sign refers to neutrinos while the minus to antineutrinos. The 
second term in Eq. (4) is CP violating. While we will use this formula to discuss our results in 
some cases, our simulation is based on the exact probabilities. 

Unless stated otherwise, we have generated our simulated data for the benchmark values in 
the first column of Table 3. These values have been chosen in conformity with the status of the 
oscillation parameters in the light of the current neutrino data [7]. The values of sin 2 26i 3 (true), 
Sqp (true) and mass hierarchy which are allowed to vary in our study, will be mentioned wherever 
applicable. For the Earth matter density, we use the PREM profile [69]. We expect to have a 
better knowledge of all the parameters mentioned in Table 3 when the Beta-beam facility comes 
up. In particular, we assume that the la error on them will be reduced to the values shown in 
the second column of Table 3 [14, 70, 71]. Therefore, we impose "priors" on these quantities, 
with the corresponding la error. The results presented in section 3 have been generated using 
the x 2 technique and numerical code described in [40, 41]. Figures in sections 4 and 5 have been 
generated using the GLoBES package [72]. For the latter simulations we do not put any priors on 
(Am 2 ! | and sin 2 2# 13 and instead add 10 year prospective disappearance data from T2K [28]. All 
other details of the \ 2 technique, as well as beam and detector specification, are taken identical in 
both numerical codes. We have made extensive checks, and the results obtained from both codes 
match to a reasonably high level of accuracy. This robustness of the results can be regarded as an 
independent cross-check within our study. 



3 Optimizing the Baseline 

We optimize the Beta-beam experiments separately with respect to the following physics outputs: 

1. The #i3 measurement reach. 

2. The mass hierarchy reach. 

3. The CP sensitivity reach. 
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Figure 2: The #13 sensitivity reach as a function of the baseline for three different values of 7 = 350 
(blue dot-dashed lines), 500 (red solid lines) and 650 (black dashed lines). The thick lines show 
the results for 8 B and 8 Li as source ions, while the thin lines give the corresponding results for 
18 Ne and 6 He. The region above the curves/within the isolated islands are permitted by the 
sensitivity criterion. 



By "reach" we refer to going to as small sin 2 2# 13 as possible. We define the performance indicators 
below, and optimize the Beta-beam experiment with respect to the baseline in this section. Note 
that we fix the 7 as well as the number of useful ion decays to their reference values in this section. 

3.1 The #13 Sensitivity/Discovery Reach 

We define two sets of performance indicators for quantifying the sensitivity of the experiment to #13. 
We call them the "#13 sensitivity reach" and the "#13 discovery reach" . The #13 sensitivity reach is 
defined as the range of sin 2 2#i3 which is incompatible with the data generated for sin 2 2#i3 (true) = 
at the 3cr CL. This performance indicator corresponds to the new sin 2 2# 13 limit if the experiment 
does not see a signal for # 13 -driven oscillations 5 . In that case, we can exclude some allowed values 
of sin 2 2# 13 , which we call our "# 13 sensitivity reach". In Fig. 2, we show the L-dependence of the 
#13 sensitivity reach of the Beta-beam experiment. The thick lines show the sensitivity for the 
8 B and 8 Li combination, while the thin lines show the corresponding sensitivities for the 18 Ne and 
6 He ions. The results are shown for three different values of 7. Since the true value of #13 is 
assumed to be zero, the data is independent of the true neutrino mass hierarchy and #cp (true), 

5 Note from Eq. (4), while the first three terms go to zero when #13 — * 0, the last term, which depends only 
on the solar parameters and #23, remains non- vanishing. Therefore, when the flux is high, i.e., for large 7 and/or 
enhanced luminosity, we expect a sizable number of events even when sin 2 29\3 (true) = 0. 
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Figure 3: The #13 discovery reach as a function of baseline for 8 B and 8 Li (left panel) and for 
18 Ne and 6 He (right panel). The pink hatched region is where a 3a discovery is possible for all 
values of Sep (true), the white boxed region is where a 3a discovery is possible for some values of 
<5cp (true), while the unshaded blank region is where its impossible to get a discovery at 3cr for 
any value of <5cp (true). The red dashed curve is the discovery reach for Sep (true) = while the 
blue dashed-dotted curve is the same for Sep (true) = tt. Here the true hierarchy is assumed to 
be normal (NH). 



but the fit depends on sgn(Am| 1 ) and Sep- We have marginalized our results over all oscillation 
parameters, including mass hierarchy and Sep- We have also marginalized over the normalization 
of the Earth matter density 6 . For the 8 B and 8 Li combination, the best #13 sensitivity is obtained 
at the magic baseline. This baseline is defined by the condition [35] 

sin(iA) ~ , (6) 

which evaluates to ^/2GFn e L(n e ) = 27r, or L ~ 7000 to 7500 km. Therefore, the second, third 
and last terms in Eq. (4) vanish at this baseline. Since the second and third terms are the CP- 
dependent terms (with the second term being CP violating), the effect of Sep is absent. Therefore, 
the correlation and degeneracies are hardly present, increasing the sensitivity of the experiment. 

The impact of the magic baseline is particularly visible for the 8 B and 8 Li combination because 
for these ions the fluxes peak at E ~ 5 — 10 GeV for 7 ~ 350 — 650. It turns out that for these 
energies, one obtains near-resonant matter effects, corresponding to A —>■ 1 in Eq. (4). Therefore, 
the flux decreases less than 1/L 2 at these energies, and is still quite substantial at the magic 
baseline. In fact, for short distances, the 1/L 2 dependence is canceled by the resonant probability 
enhancement, as can be read off Eq. (4). For the 18 Ne and 6 He combination, the fluxes peak 
around E ~ 1.0 — 2.5 GeV for 7 ~ 350 — 650. For such low energies, matter effects are small, even 
for very long baselines. More importantly, for E ~ 1.0 — 2.5 GeV, the oscillatory factor peaks 
at L ~ 500 — 1250 km if one assumes Am^ = 2.5 x 10~ 3 eV 2 . Therefore, for this ion pair, the 

6 In all results given in this paper, we have done full marginalization over hierarchy, all oscillation parameters 
and the normalization factor of the Earth matter density distribution. 
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minimum in the sin 2 2#i 3 sensitivity comes at the baseline where we expect the first oscillation 
maximum. For 7 = 650, 500 and 350, the best sensitivity comes at L — 1250 km, 890 km and 
680 km, respectively. Beyond this baseline, both the flux and the probability fall, resulting in a 
sharp loss of the sin 2 2#i 3 sensitivity. Note the isolated regions for 18 Ne and 6 He and 7 = 500 and 
650, which are also incompatible with sin 2 2^13 (true) = 0. The gap is mainly an artifact of the 
presence of clone solutions at these smaller baselines and it might be breached by the combination 
of a higher flux, better energy resolution, etc., leading to a much better sensitivity. 

In Fig. 3, we show the "discovery reach" for sin 2 2# 13 (true). This performance indicator is 
defined as the range of sin 2 2# 13 (true) values which allow us to rule out sin 2 2# 13 = at the 
3o" CL. Since the data are now generated for a non-zero sin 2 2^13 (true), there is a <5cp (true) 
dependence and a true mass hierarchy dependence. The discovery reach for <5cp (true) = is 
shown by the dashed curve in Fig. 3, and the discovery reach for <5cp (true) = 7r is shown by the 
dashed-dotted curve. For each <5cp (true), one obtains a corresponding such curve. To show the 
impact of this <5cp (true) dependence of the discovery potential, and to illustrate explicitly the 
increase in the "risk factor" coming from our lack of knowledge of <5cp (true), we present in Fig. 3 
a band marked by boxes showing the entire range of sin 2 2#i 3 (true) values corresponding to all 
possible values of <5cp (true). Fig. 3 has been drawn for the true normal hierarchy. 

The way to interpret this figure is as follows: At a given L, one will discover sin 2 2# 13 for any 
5cp (true) at the upper limit and beyond (the pink hatched region), whereas there is no value of 
Sep (true) for which one can discover sin 2 2#i 3 below the lower limit (the unshaded region). Within 
the band, the fraction of Sep (true) , which allows for a discovery, increases as one approaches the 
upper limit. Therefore, the upper edge of this band gives the most conservative discovery reach. 
We will take this as our final discovery reach at a given baseline. The lower edge of the band 
shows the best possible case. The left panel of Fig. 3 is computed for 8 B and 8 Li, while the right 
panel is computed for 18 Ne and 6 He . All results in this figure have been computed for 7 = 500. 
Again we note that, at the magic baseline, our results are CP independent, and the band reduces 
to a point since the discovery reach is independent of 5 C p (true) at this baseline 7 . For the 8 B and 
8 Li combination the best discovery reach comes at L — 7600, which is the magic baseline. It 
is noteworthy though that for these ions the best possible case, given by the lower edge of the 
band, comes at a lower baseline of L ~ 700 km. For the 18 Ne and 6 He pair, since matter effects 
are low, there is only a local minimum at the magic baseline. In fact, the best performance of 
the experiment comes at L = 900 km, which is approximately the baseline where we have the 
oscillation maximum for 7 = 500. Again we note the appearance of islands inside the band for 
the 18 Ne and 6 He ions. For regions inside these islands, the discovery of sin 2 26>i 3 is independent 
of Sep (true). 

We have also computed the sin 2 2#i 3 (true) discovery reach for the true inverted hierarchy. 
Since the figures look very similar to Fig. 3, we do not show them explicitely. However, note 
that in the inverted case, the curves for 5qp (true) = and tt interchange their roles. It is not 
surprising that the inverted hierarchy performs similar to the normal one for the beta beams, since 
the neutrino and antineutrino event rates are very similar. Our chosen number of useful ion decay 
is about a factor of three higher for antineutrinos than neutrinos, which is compensated by the 

7 For the 18 Ne and 6 He combination, there is a small width even at the magic baseline. In this case sin 2 29is (true) 
is large, and Eq. (4) has to be expanded to higher order. The magic baseline condition Eq. (6) may not hold for 
higher order terms. This was also noted and pointed out in [41]. 
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Figure 4: sin 2 2#i3 (true) reach for determining sgn(Am| 1 ) at 3a CL as a function of the baseline 
with 8 B and 8 Li (left panel) and 18 Ne and 6 He (right panel) taken as the source ions. The pink 
hatched region is where the wrong hierarchy can be ruled out for all values of <5cp (true), the 
white boxed region is where the hierarchy determination is possible for some values of 5qp (true), 
while the unshaded blank region is where its impossible to determine the hierarchy for any value 
of <5cp (true). The red dashed (blue dashed-dotted) curves show the hierarchy sensitivity for 
5qp (true) = (Sep (true) = ir). The true hierarchy is assumed to be normal (NH). 



higher neutrino cross sections. 

3.2 The sgn(Ani3 1 ) Sensitivity Reach 

We define the mass hierarchy sensitivity as the range of sin 2 2# 13 (true) for which the wrong 
hierarchy can be excluded at the 3a CL. We show our results as a function of the baseline in 
Fig. 4. As before, the left panel is for the 8 B and 8 Li case, and the right panel for 18 Ne and 6 He. 
In addition, we show the risk with respect to <5cp (true) as a band marked by boxes, where the 
lower edge corresponds to the best possible reach (obtainable for only some specific Sep (true)), 
and the upper edge to the conservative case (valid irrespective of the value of Sqp (true)). That 
means that the hierarchy can be determined for any sin 2 2# 13 (true) above the upper end of the 
band. In all panels, we also show the curves corresponding to 6qp (true) = and 6qp (true) = tt, 
for illustration. For the 8 B and 8 Li case, we find that the best sensitivity to the mass hierarchy 
comes at the magic baseline. The reason for this is basically the same as in the previous subsection: 
The near-resonant matter effects lead to a large number of events, and the resonant behavior is 
only present for one hierarchy (normal or inverted). It is, therefore, possible to have very large 
matter dependent oscillations for L ^ 4000 km. However, the effect of Sqp (true) could wash 
away the sensitivity to sgn(Am| 1 ). For example, for L ~ 4000 km, we do not obtain a very 
good sensitivity for even the best case. If one takes into account all possible <5 CP (true) values, the 
sensitivity becomes deteriorated significantly. At the magic baseline, the dependence on Sqp (true) 
is reduced. Therefore, this baselines provides the best choice to determine the mass hierarchy. 
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The right panel in Fig. 4 corresponds to the 18 Ne and 6 He case. The sensitivity to the mass 
hierarchy is rather poor for the values of 7 we have adopted here because of the low energies off 
the matter resonance (the energies are about a factor of 3.5 lower than for the 8 B and 8 Li pair). 
Since the mass hierarchy determination crucially depends on matter effects, we have very poor 
sensitivity for this performance indicator, at least for the values of 7 considered here. We will see 
in the next section that this set of ions could start giving comparable sensitivity only when the 7 
is increased by a factor of three. For 7 = 500, the best sensitivity comes at L ~ 2000 km. 

Again, we have tested the true inverted hierarchy case, and we have not found any significant 
qualitative or quantitative differences. 

3.3 The CP Sensitivity Reach 

We next discuss the reach of the experiment to CP violation as a function of L. We define the 
sensitivity to (maximal) CP violation as the range of sin 2 2# 13 for which CP conservation (i.e., 
Sep = and ir) can be excluded at the 3a confidence level irrespective of the (fit) hierarchy. The 
results are shown in Fig. 5, where we generate the data either for <5cp (true) = n/2 (left panels) 
or 37r/2 (right panels). The upper panels are computed for the 8 B and 8 Li pair, while the lower 
panels are for the 18 Ne and 6 He pair. We show results for three choices of 7. Let us focus on 
Sep (true) = 7r/2 (left panels) first. The best choice for the baseline is around L = 500 — 1500 
km, depending on the choice of ions and 7. For example, for the 8 B and 8 Li combination, the 
optimum is found at around the FNAL-Homestake baseline L = 1290 km. However, the absolute 
performance for 8 B and 8 Li is worse than for the 18 Ne and 6 He combination. For 37r/2 (right 
panels), the sensitivity becomes worse in both cases due to the impact of degeneracies. Ignoring 
the gaps, there is still a substantial CP violation reach for 18 Ne and 6 He for both Sep (true) = 7r/2 
and 3n/2. For example, if Double Chooz constrains sin 2 2^13 to values smaller than about 0.04, 
a short baseline L ~ 300 km together with 7 > 500 might be sufficient for the CP violation 
measurement because the gap can be excluded. The optimal baseline for maximum sin 2 2#i 3 
reach, on the other hand, ranges between L ~ 600 for 7 = 350 to L = 1000 km for 7 = 650. For 
the inverted hierarchy, we obtain very similar figures to Fig. 5, but the role of 5qp (true) = %/2 
and 5qp (true) = 37r/2 is exchanged. 

3.4 Comparison with a Megaton Water Detector Set-up 

In this paper, we have restricted ourselves to magnetized iron calorimeters as the far detector 
technology. Detector characteristics such as size, efficiency, energy threshold, energy resolution 
and backgrounds affect crucially the choice of the source ions and their corresponding boosts. 
Energy threshold is expected to be lower for water Cerenkov detectors, totally active scintillator 
detectors, and liquid argon TPC. As an additional complication, the reconstructed energies of the 
background events can be different. Therefore, the optimization of the Beta-beam depends on the 
detector technology as well. In this paper, we have restricted ourselves to a single detector option 
for all cases, in order to outline the optimization in terms of boost factor, baseline, luminosity, 
and especially the ion pair. A comparison of different detector technologies will follow [74] . Here 
we compare our results to a water Cerenkov detector set-up for the sake of illustration. We 
use as a benchmark the set-up from Ref. [46], in which 18 Ne and 6 He ions with 7 = 350, L = 
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Figure 5: sin 2 29 13 (true) reach for sensitivity to maximal CP violation as a function of L, for three 
different values of 7. The upper panels are for 8 B and 8 Li taken as the source ions and the lower 
panels are for 18 Ne and 6 He as the source ions. The left panels are for Sep (true) = 7r/2, and the 
right panels for 5qp (true) = 37r/2. In all true normal hierarchy has been assumed. 
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(L^Magic, 8 B+ 8 Li) 


(L^Magic, 8 B+ 8 Li) 


(L^Magic, 8 B+ 8 Li) 


(L=700 km, 18 Nc+ 6 He) 


Benchmark 

water Cerenkov 
7 = 350 
L = 730 km 
18 Ne+ 6 He 


5.7 x 10~ 4 


8 x 10~ 4 


1.6 x 10~ 2 


5.2 x 10~ 5 



Table 4: Comparison of the conservative sensitivity reaches in sin 2 2# 13 for the different perfor- 
mance indicators. Conservative is defined as the reach in sin 2 2#i3 irrespective of 8cp (true) and 
true mass hierarchy. For CP violation we give the reach for Sep (true) = n/2. The upper row 
gives the sensitivity reaches for the optimal set-ups identified in this paper, which use ICAL-type 
magnetized iron detectors. The value of L and the ion source are shown in parentheses. Lower 
row gives the sensitivity for the benchmark set-up with water Cerenkov detector. The exposures 
correspond to 50 kton of iron and 500 kton of water with five years of neutrino and antineutrino 
runs each. 

730 km, and water Cerenkov detector with 4.4 Mt x y statistics were used as one option. Such 
a detector could be placed at the Gran Sasso laboratory in Italy or at Canfranc in Spain. If 
the source is located at CERN the baselines are 730 km and 650 km, respectively. In order to 
compare with the results of our optimal set-ups, we have re-computed the sensitivities of this 
benchmark set-up with GLoBES [72] for the same input assumptions as in this study and a total 
luminosity corresponding to 500 kton of water with five years of neutrino and antineutrino runs 
each 8 . This set-up returns a sin 2 2^13 sensitivity of 5.7 x 10~ 4 , while the #13 discovery reach 
ranges between 4.8 x 10~ 5 and 8.0 x 10 -4 depending on Sqp (true). For two specific choices of 
Sep (true) = and 7r the discovery reaches are 3.1 x 10 -4 and 8.0 x 10~ 4 respectively. Normal mass 
hierarchy discovery reach ranges between 2.4 x 10~ 3 and 1.6 x 10~ 2 depending on 5 C p (true). For 
5cp (true) = and 71 the hierarchy sensitivity reaches are 1.0 x 10~ 2 and 4.0 x 10~ 3 respectively. 
Maximal CP violation can be established if sin 2 2^13 (true) > 5.2 x 10~ 5 for <5cp (true) = n/2 
and sin 2 2# 13 (true) > 5.5 x 10~ 5 for Sqp (true) = 37r/2. All numbers at the 3a confidence level, 
and we have not found any disconnected regions for these performance indicators. For a direct 
comparison with the optimal set-ups identified in this paper, we present in Table 4 the "most 

8 Especially, the different values of 6*12 and Am^ compared to Ref. [46] do have some impact on the sensitivities. 
The experiment description is taken from the current GLoBES distribution, updated with the luminosity numbers 
for this study. 



17 



conservative" sensitivity reaches in sin 2 29\s (true) for the four types of performance indicators we 
have defined. "Most Conservative" is defined as the reach after allowing for all possible 5 C p (true) 
and true mass hierarchy This corresponds to a "no-risk" situation. The upper row gives the 
sensitivity reaches for our optimized set-ups for the different performance indicators. The set- 
up concerned is defined within the parentheses. The lower row gives the corresponding reaches 
for the set-up involving water Cerenkov detector defined above. One can see that the water 
Cerenkov set-up comprehensively outperforms our optimal set-up in CP violation, while mass 
hierarchy measurement is clearly done better with 8 B and 8 Li ions at the magic baseline. The 
#13 measurement is better in water Cerenkov set-up, though our optimal set-up is only slightly 
worse. Also, it is unclear if the 50 kton iron calorimeter is comparable to the 500 kton water 
Cerenkov detector (fiducial volume) in terms of cost. 

4 Optimizing the Baseline and Gamma 

In this section, we will allow for larger values of 7 and show results in the 7 — L plane, while we 
still fix the number of useful ion decays to its reference value. In section 2, we have compared the 
different ion pairs for the same fixed 7. This however, implies that the peak energies, and therefore 
the physics (matter effects etc.) is different. Here we wish to show the projected sensitivity for 
the two sets of isotopes when they produce beams with similar peak energies, where we leave both 
L and 7 as free parameters. The motivation is to have similar effect of the oscillations and matter 
effects for both the sets of ions. The end-point energies of 18 Ne and 6 He are about 3.5 times 
smaller than those of 8 B and 8 Li . Therefore, with 7 for 18 Ne and 6 He scaled up by a factor of 
roughly 3, we expect almost the same oscillated spectral shape for the two cases. Hence we show 
the results for 7 < 1000 for 8 B and 8 Li, and for 7 < 3000 for the 18 Ne and 6 He. We stress that we 
allow for these prohibitively large values of 7 for 18 Ne and 6 He in order to compare set-ups with 
same physics output. While they might appear to be unrealistic, they can in principle be achieved 
by an accelerator as large as the LHC [73]. The Tevatron might also be used to give large boost 
factors. 

4.1 The 0i3 Sensitivity 

In order to discuss the #13 performance, we use the sin 2 2#i3 sensitivity reach in this section. The 
sin 2 2#i3 sensitivity represents the largest fit sin 2 26>i 3 which fits a sin 2 2#i 3 (true) = 0. Therefore 
in this section, we do not take into account any disjoint regions, such as the ones which appeared 
in Fig. 2. We reiterate that the sin 2 2# 13 sensitivity is independent of the true mass hierarchy 
and Sep (true). We show in Fig. 6 the 3<r sin 2 2# 13 sensitivity reach as a function of baseline 
L and Lorentz factor 7. The curves in this figure represent the 3a C.L. contours for different 
values of sin 2 29 13 . The contours are spaced by 0.2 in log 10 (sin 2 29 13 ), and the numbers are shown 
in the figure for some of them. The diamonds mark the absolute optimum within each plot; the 
corresponding values are given in the figure caption. The upper row corresponds to 8 B and 8 Li and 
the lower row to 18 Ne and 6 He. The left column shows the systematics limit only, where we have 
kept all parameters including the mass hierarchy and <5cp = fixed in the fit. The right column 
shows the final sensitivity limit after including correlations and degeneracies. From Section 2.1, 
we have learned that a gamma range for 18 Ne and 6 He about three times larger than for 8 B and 
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Figure 6: The 3a sin 2 2#i3 sensitivity as a function of baseline L and boost factor 7. The sin 2 26*13 
sensitivity represents the largest possible sin 2 26*i 3 , which fits data simulated at sin 2 26*i 3 (true) = 0. 
The upper row corresponds to 8 B and 8 Li, the lower row to 18 Ne and 6 He. The left column shows 
the systematics limit only (i.e., the oscillation parameters are fixed with <5cp = 0), whereas the 
right column shows the final limit including correlations and degeneracies. The contours are spaced 
by 0.2 in log 10 (sin 2 2^13), where the numbers are given for some of these. The diamonds mark the 
absolute optimum within each plot, which are 10~ 3 ' 66 , 10 -3 00 , 10 -4 ' 66 , and 10 -4 ' 21 , respectively, 
from the upper left to the lower right. On the right axes of the plots, an energy scale is attached 
which corresponds to the mean peak energy jEq with E being the mean of the endpoint energies 
for the isotope pair. 
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8 Li gives similar neutrino peak energies, which is illustrated by the right vertical axes in the plots. 
This feature is reflected in our choice of the 7 ranges in Fig. 6. Note that in both the upper and 
lower rows of Fig. 6 we have used the same isotope decay rates, the same running times, and the 
same detector simulation. 

Let us first focus on the systematics limit shown in the left column. The two panels for the 
different isotopes look qualitatively very similar, but the absolute performance is much better for 
18 Ne and 6 He for 7 values three times larger. This is because while we have tuned the oscillation 
probability for both isotopes to be the same, the flux increases as 7 2 , and therefore the event 
rate for the 18 Ne and 6 He combination scales up by a factor of about nine for the three times 
higher boost factors. According to Eq. (2), the absolute sensitivities for the two isotopes would 
be more or less identical if the 7 for 18 Ne and 6 He was turned up by a factor of about 3.5 and 
the source luminosity for 8 B and 8 Li was enhanced by a factor of about 12. This means that 
primary difference between the upper and lower rows in Fig. 6 is the different event rate, but not 
the spectral shape and energies. We can also note from the figure that for a give set of source 
ions and for a given L and Np, the sensitivity of the experiment increases as 7 is increased. For 
the best sensitivity reach, the baseline has then to be adjusted accordingly. The figure shows that 
the best choice for L roughly corresponds to tuning L/7 ~ 2.6 for 8 B and 8 Li and L/7 ~ 0.8 for 
18 Ne and 6 He. 9 We show the lines for these conditions in Fig. 6 in the upper left panel and lower 
panels 10 . Note that the slopes of these lines are different by about a factor of / ~ 3.5 
because the boost factors are related by Eq. (2) in order to obtain the same neutrino energies. 
The optimal regions appear slightly to the left of these lines in the figure because statistics are 
higher at lower L. 

Once the correlations and degeneracies are included (c./., right column of Fig. 6), the optimum 
baseline changes qualitatively for 8 B and 8 Li, but not for 18 Ne and 6 He. One can read off from 
the upper right panel that the magic baseline becomes the optimum baseline for about 7 > 350, 
whereas for 7 < 350 the shorter baseline is preferred. The reason for this was also discussed in 
the previous section. For 7 > 350 one gets a peak neutrino energy greater than 5 GeV. For these 
energies one can obtain very large matter effects for the very long baselines and hence the oscillation 
probability becomes larger here. Therefore, the event rate in this regime is improved due to a 
combination of large matter driven oscillation probability as well as increased flux driven by larger 
7. Most importantly, close to the magic baseline, the effect of <5cp is almost negligible. Therefore, 
once the correlations and degeneracies are taken into account, this becomes the deciding factor 
which ensures that for 7 > 350, the magic baseline emerges as the most optimal baseline choice 
for the #13 sensitivity. For 7 < 350, both the oscillation probability as well as the flux collimation 
are small, and therefore the longer baselines suffer due to the 1/L 2 suppression factor for the event 
rates. As a result, despite being free from problems of correlations and degeneracies, the magic 
baseline looses sensitivity compared to the shorter baselines. Also note that the sensitivity for 
7 < 350 is rather poor, and therefore it is desirable to do the experiment at 7 > 350. 

The situation is very different for the 18 Ne and 6 He case, where the shorter baseline (which 
roughly satisfies L/7 ~ 0.8) is always the better choice for the sin 2 2^13 sensitivity for any value 
of 7. Note that we have adjusted 7 to keep the same peak energy (shown by the dashed lines) 

9 In fact, the unit of L/7 is km, which we do not put explicitly. 

10 Sincc for 8 B and 8 Li we do not get the optimum at the oscillation maxima obeying L/7 ~ 2.6 once correlations 
and degeneracies are taken, we do not show the line in the upper right panel. 
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and hence we have the same spectral shape and the same effect of the oscillation probability as 
compared to 8 B and 8 Li. However, the 7 here is about three times larger, and that becomes the 
overwhelming deciding factor for the most suitable baseline. The enhancement in the flux due to 
the beam collimation in this case ensures that the shorter baselines have a high enough statistics to 
handle the problem of parameter correlations and degeneracies. Therefore, the sin 2 2#i3 sensitivity 
is better than that of the magic baseline, which is free from a <5cp dependence but suffers from the 
1/L 2 suppression. Therefore, the best sensitivity roughly follows the L/7 = 0.8 line. We remind 
the reader that the isolated regions obtained for 7 = 500 and 650 in Fig. 2 also correspond to 
L/7 ~ 0.8. For small values of sin 2 2# 13 , that are relevant here, the Sqp dependent terms in the 
fit are extremely important at values of L close to L oscmax /2. This corresponds to L/7 ~ 0.8 for 
18 Ne and 6 He. 

In the following section, we will compare the long (magic) with the short baseline. For the 
long (magic) baseline, we will choose L = 7500km (optimum). For the short baseline we will take 
the local optimum obtained, which obeys L/7 ~ 2.6 for 8 B and 8 Li and L/7 ~ 0.8 for 18 Ne and 
6 He. 

4.2 The sgn(Am3 1 ) Sensitivity Reach 

In Fig. 7, we show the sin 2 2#i3 (true) reach for the sensitivity to the normal mass hierarchy (3a) 
as a function of baseline L and boost factor 7. The sgn(Am| 1 ) reach represents the minimum 
sin 2 2^13 (true) above which the mass hierarchy will be discovered for any sin 2 26*i 3 (true) (i.e., 
there are no gaps in the sensitivity). The upper row corresponds to 8 B and 8 Li, the lower row to 
18 Ne and 6 He. The left column is computed for 5 CP (true) = tt/2, whereas the right column is 
for Sep (true) = 3tt/2. Note that the sin S term in Eq. (4) is positive for neutrinos and negative 
for antineutrinos. This means that the simulated neutrino rate is larger for 5qp (true) = tt/2 and 
smaller for Sqp (true) = 3n/2. Therefore, statistics is better for <5cp (true) = tt/2, and we expect 
a better mass hierarchy sensitivity. It is for this reason we have chosen to illustrate the sgn(Am§ 1 ) 
sensitivity reach for 5qp (true) = n/2 and 3ir/2 only in Fig. 7, which represent two cases close to 
the best case and worst case. These two cases would change their role if one used a true inverted 
hierarchy instead, i.e., the best performance would be close to Sep (true) = 3n/2, and the worst 
close to 5qp (true) — tt/2. 

For both isotope pairs, there are two main observations for the mass hierarchy measurement: 

1. Longer baselines are preferred with the optimal sensitivity reach appearing at L close to the 
magic baseline. This is not surprising since for longer baselines the matter effect contribution 
becomes larger, which allows to discriminate between normal and inverted hierarchy. The 
sensitivity gets better at the magic baseline since the probability at the magic baseline is 
free of 5qp related correlations and degeneracies. 

2. Higher boost factors are preferred (at least within the shown ranges). First of all, the event 
rate increases as 7 increases. However, matter effects also increase as one approaches the 
mantle resonance energy at about 7 GeV. Therefore, sufficiently high energies are needed to 
observe the mass hierarchy discriminating matter effects. 

In summary, we find that the optimal choice for determining the hierarchy is L ~ 7000 — 
9000 km, and 7 > 200 for 8 B and 8 Li or 7 > 750 for 18 Ne and 6 He. Of course, for the same 
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Figure 7: The sin 2 2#i3 (true) reach for the sensitivity to the normal mass hierarchy (3a) as a 
function of the baseline L and boost factor 7. The sin 2 2# 13 (true) reach represents the minimum 
sin 2 26*13 (true) above which the mass hierarchy will be discovered for any sin 2 26 ) 13 (true) (i.e., 
there are no gaps in the sensitivity). The upper row corresponds to 8 B and 8 Li, while the lower 
row to 18 Ne and 6 He. The left column is computed for Sqp (true) = tt/2, whereas the right column 
is for 5qp (true) = 37r/2. The contours are spaced by 0.2 in log 10 (sin 2 26>i 3 (true)), where the 
numbers are given for some of these. The diamonds mark the absolute optimum within each plot, 
which are 10 -3 ' 46 , 10 -3 ' 36 , 10~ 4 ' 29 , and 10~ 4 ' 10 , respectively, from the upper left to the lower right. 
On the right axes of the plots, an energy scale is attached which corresponds to the mean peak 
energy 7^0 with Eq being the mean of the endpoint energies for the isotope pair. 
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energy, the absolute performance is better for 18 Ne and 6 He as long as one can create boost factors 
which are three times larger. 

4.3 The CP Sensitivity Reach 

In order to discuss the CP sensitivity, we show in Fig. 8 the sin 2 2#i 3 (true) reach for the sensitivity 
to maximal CP violation (3a) as a function of the baseline L and boost factor 7. The data are 
generated for normal hierarchy and <5cp (true) = 7r/2 (left panels) or 3n/2 (right panels). We show 
in the figures the minimum sin 2 2# 13 (true) above which maximal CP violation will be discovered 
for any sin 2 2^13 (true). In these figures we discard the lower allowed islands in Fig. 5 (lower right 
panel) and consider only the upper regions. The upper row corresponds to 8 B and 8 Li, the lower 
row to 18 Ne and 6 He. Obviously, for the CP violation sensitivity, a shorter baseline is a must, 
because the magic baseline is not sensitive to Sep- Hence, we only show baselines up to 5000 
km in this figure, since there is no sensitivity for longer baselines. As for the mass hierarchy 
reach, here too the performance for the normal hierarchy is best for Sep (true) = tt/2 and worst 
for Sep (true) = 37r/2. However, their roles change for the inverted hierarchy. Note that the 
poor sensitivity close to Sep (true) ~ 37r/2 mainly comes from unresolved degeneracies due to 
poor statistics. Combining data of this set-up with a second (much longer) baseline could help 
in resolving the degeneracies and improving the sensitivity. Here we focus on the left column of 
Fig. 8 for the following discussion. 

A comparison of the upper left and lower left panels of Fig. 8 reveals qualitative differences 
between the sensitivities coming from the two pairs of isotopes. We observe that the 18 Ne and 
6 He combination is far superior for probing the CP phase. We had already seen this feature in 
the earlier section when we had compared the CP sensitivity reach of the two sets of isotopes for 
the same 7. In fact, the sensitivity to maximal CP violation improves with 7 for 18 Ne and 6 He, 
and we find that the best case shown by the diamond appears for the highest 7 we have taken. 
For 8 B and 8 Li the sensitivity in general is comparatively poorer and does not scale with 7. In 
fact, the 7 dependence of the sensitivity is rather weak with the best CP sensitivity coming for 
7 ~ 650. As far as the possible baselines are concerned, for <5cp (true) = n/2 the optimal baselines 
roughly follow the L/7 = 2.6 line for 8 B and 8 Li, and the L/7 = 0.8 line for 18 Ne and 6 He. 

5 The Impact of Luminosity 

In this section, we study the impact of increasing the overall number of events through either 
increasing the size of the detector, the exposure time, detector efficiency, or the number of useful 
ion decays. In order to discuss this, we introduce a luminosity scaling factor multiplying the 
overall luminosity (useful isotope decays x running time x detector mass x detector efficiency) 
for both neutrinos and antineutrinos. Note that the luminosity scaling factor corresponds to a 
reference luminosity, i.e., 1.1 x 10 18 (year _1 ) x5(year) x50(kton) x0.76 for the neutrino beam. 
For the antineutrino beam we use 2.9 x 10 18 useful decays per year. We attempt to determine the 
optimal value for the luminosity, 7 and L. We fix 7 at certain benchmark values and compare 
the performance of short baselines given by L/7 ~ 2.6 for 8 B and 8 Li and L/7 ~ 0.8 for 18 Ne and 
6 He, with the magic baseline. For comparison between the isotopes, we use both the approaches 
discussed in section 2.1. That is, we compare the physics reach of 8 B and 8 Li with that of 18 Ne and 
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Figure 8: The sin 2 2#i3 (true) reach for the sensitivity to maximal CP violation (3a) as a function 
of the baseline L and boost factor 7. The sin 2 2#i3 reach represents the minimum sin 2 2#i3 (true) 
above which maximal CP violation will be discovered for any sin 2 2# 13 (true) (i. e., there are no gaps 
in the sensitivity). The upper row corresponds to 8 B and 8 Li, the lower row to 18 Ne and 6 He. The 
left column is computed for 5qp (true) = tt/2, whereas the right column is for 5qp (true) = 37r/2. 
The contours are spaced by 0.2 in log 10 (sin 2 2#i3 (true)), where the numbers are given for some of 
these. The diamonds mark the absolute optimum within each plot. Here a true normal hierarchy 
is assumed. The diamonds mark the absolute optimum within each plot, which are 10~ 2 ' 82 , 10~ 2 ' 14 , 
10 -4 ' 50 , and 10 -2 ' 79 , respectively, from the upper left to the lower right. On the right axes of the 
plots, an energy scale is attached which corresponds to the mean peak energy 7^0 with Eq being 
the mean of the endpoint energies for the isotope pair. 
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6 He, both at same fixed values of 7 corresponding to the same input, as well as with 7 for 18 Ne and 
6 He scaled by a factor of about 3.5 to get the same peak neutrino energy, and hence the same 
neutrino energies. We use the same definition of the performance indicators as in Section 4, i.e., 
we demand that there is sensitivity for all sin 2 26 13 larger than the given sensitivity limits (which 
excludes the regions separated by the gaps). 

5.1 The #13 Sensitivity 

In Fig. 9, we show the sin 2 2#i3 sensitivity (3<r) as a function of the luminosity scaling factor. The 
panels represent the different isotopes and different 7 values as given in the captions. The dot- 
dashed curves correspond to the magic baseline (MB) with L = 7500 km fixed, the solid curves to a 
short baseline with an L/7 depending on the isotope. The upper row is for 8 B and 8 Li. The middle 
and lower rows are for 18 Ne and 6 He, with the middle row for same 7 as 8 B and 8 Li, whereas the 
lower row is for 7's scaled up by a factor of about 3.5. That means that the middle row represents 
the same accelerator effort as the upper row in terms of 7, and the lower row represents similar 
neutrino energies to the upper row. 

There are a number of interesting observations from Fig. 9. First, for the shorter baseline, 
statistics are crucial for resolving the degeneracies, and we obtain a sudden enhancement of the 
sensitivity at some value of the luminosity scaling factor. This point is visible as the edge in 
the luminosity scaling, where the degenerate solution is ruled out at the 3<r C.L. For the magic 
baseline, degeneracies are hardly relevant and we find that the sensitivities exhibit an almost power 
law scaling with statistics. Second, the curves for the shorter baseline cross the ones for the magic 
baseline twice in almost all the panels in the upper and lower rows. For 8 B and 8 Li, the standard 
assumed luminosity (luminosity scaling factor one) is typically in the window where the magic 
baseline performs better. Only for the 7 = 350 case, the shorter baseline is comparable with the 
magic baseline for the standard luminosity. For 18 Ne and 6 He, the short baseline is typically better 
for the standard luminosity and smaller 7 (middle row). However, for 18 Ne + 6 He and 7 3> 1000, 
already a factor of two loss in luminosity makes the magic baseline the better choice. 

Let us now come back to the conditions in Eq. (1). Since we have chosen the upper and lower 
rows in Fig. 9 such that the gammas scale inverse to the endpoint energies, this formula indicated 
that the same physics should be obtained for about a factor of 12 difference in luminosity. Let us 
pick a simple feature where we could test this conclusion: take a look at the edge of the sensitivity 
jumps for the short baselines. In the lower row, this sensitivity jumps happen for about one order 
of magnitude less luminosity than in the upper row, which confirms the expectation. This means 
that a factor of 12 in the number of useful ion decays is indeed needed to reproduce the same 
physics if one uses 8 B and 8 Li instead of 18 Ne and 6 He. However, this comes for a price: A factor 
of 3.5 higher gamma is needed for 18 Ne + 6 He than for 8 B + 8 Li. 

5.2 The sgn(Amg 1 ) Sensitivity Reach 

In Fig. 10, we show the sin 2 2# 13 reach for the sensitivity to the normal mass hierarchy (3er) as a 
function of the luminosity scaling factor. The panels represent the different isotopes and different 
7 values as given in the captions. The dashed-dotted curves correspond to the long magic baseline 
"MB" with L = 7500 km fixed and very little dependence on the true 5cp, the other two sets of 
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Figure 9: The sin 2 2t9 13 sensitivity (3a) as a function of a luminosity scaling factor. The luminosity 
scaling factor multiplies the overall luminosities (useful isotope decays x running time x detector 
mass x detector efficiency) for both neutrinos and antineutrinos. The panels represent the different 
isotopes and different 7 as indicated in the captions. The green dashed-dotted curves correspond 
to the magic baseline "MB" with L = 7500 km fixed, the red solid curves to a short baseline with 
an L/7 depending on the isotope (c.f., lines for fixed L/7 in Fig. 6). A true normal hierarchy is 
assumed. 
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Figure 10: The sin 2 2#i3 reach for the sensitivity to the normal mass hierarchy (3a) as a function 
of a luminosity scaling factor. The luminosity scaling factor multiplies the overall luminosities 
(useful isotope decays x running time x detector mass x detector efficiency) for both neutrinos 
and antineutrinos. The panels represent the different isotopes and different 7 as indicated in 
the captions. The green dashed-dotted curves correspond to the long magic baseline "MB" with 
L = 7 500 km fixed and very little dependence on the true Sep, the other two sets of curves to a 
short baseline with an L/7 depending on the isotope (c.f., lines for fixed L/7 in Fig. 6). The red 
solid curves are computed for a true <5cp = vr/2, the blue dashed curves for a true 5qp = 3ir/2. 
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curves to a short baseline with an L/7 depending on the isotope (c.f., lines for fixed L/7 in Fig. 6). 
The solid curves are computed for a 5 C p (true) = 7r/2, the dashed curves for a 5 C p (true) = 3n/2. 

For the mass hierarchy sensitivity, the magic baseline exhibits almost a power law scaling with 
statistics, with the best absolute performance in all cases for all luminosities in the upper and 
lower rows. For the shorter baselines, the performance depends crucially on the true 5qp- For 
Sep — 7r/2, the scaling behaves similar to the magic baseline, and there are only very few jumps. 
That means that changes in the chosen standard luminosity (scaling factor one) do not have a 
strong effect. Only for the 18 Ne and 6 He combination with very high 7, increasing the luminosity 
resolves the (wrong hierarchy) intrinsic degeneracy and this improves the performance of the set- 
up such that the sensitivities approach the ones obtained for the magic baseline. For 5qp — 37r/2, 
the mass hierarchy sensitivity is basically not present for 8 B and 8 Li, and for 18 Ne and 6 He for 
7 <C 2000. In summary, for the mass hierarchy sensitivity, the magic baseline is a safe choice, 
independent of the luminosity. Similar sensitivity can be achieved by the shorter baseline only for 
extreme choices of 7 and Np. 

5.3 The CP Sensitivity Reach 

The sin 2 2# 13 reach for the sensitivity to maximal CP violation (3a) as a function of the luminosity 
scaling factor is shown in Fig. 11. The luminosity scaling factor multiplies the overall luminosities 
(useful isotope decays x running time x detector mass x detector efficiency) for both neutrinos 
and antineutrinos. The panels represent the different isotopes and different 7 values as given in 
the captions. For all curves, the short baselines are used, i.e., L/7 = 2.6 for the upper row, 
and L/7 = 0.8 for the middle and lower rows, because there is no CP violation sensitivity at 
the magic baseline. The solid curves are computed for 5 CP (true) = tt/2, the dashed curves for 
5qp (true) = 3n/2. 

The scalings are very similar to the mass hierarchy sensitivity: For 5qp — tt/2 there are no 
jumps, and for <5cp — 3n/2 the intrinsic degeneracy can be resolved above a certain threshold 
luminosity (because it is then lifted over the x 2 — 9 line). All curves look qualitatively very 
similar, no matter what 7 or isotope is used. However, The absolute performance is better for 
18 Ne + 6 He even for the same 7. 

As far as the chosen reference luminosity (scaling factor one) is concerned, the 8 B and 8 Li beam 
is operated in a region where substantial luminosity changes will not affect the result. However, 
for 18 Ne and 6 He, already a factor of two to three luminosity increase would be sufficient to resolve 
the intrinsic degeneracy and to boost the performance at 5qp (true) = 3tt/2 for the lower row. 
Note that such a boost could also be achieved by a synergistic degeneracy resolver, such as the 
magic baseline as second baseline. 

5.4 Comparison with a Neutrino Factory 

As far as the comparison to a Neutrino Factory is concerned, our Figs. 6 (#13), 7 (sgn(Am| 1 )), 
and 8 (maximal CP violation) correspond to Figs. 3 (6*13), 6 (sgn(Am| 1 )), and 5 (maximal CP 
violation) in Ref. [38] for a Neutrino Factory, which means that they allow for a direct qualitative 
comparison. While for the beta beams, the boost factor is on the vertical axis, for a Neutrino 
Factory, the muon energy is on the vertical axis. Both are proportional to the peak neutrino 
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Figure 11: The sin 2 29i3 reach for the sensitivity to maximal CP violation (3a) as a function of a 
luminosity scaling factor. The luminosity scaling factor multiplies the overall luminosities (useful 
isotope decays x running time x detector mass x detector efficiency) for both neutrinos and 
antineutrinos. The panels represent the different isotopes and different 7 as given in the captions. 
For all curves, the short baselines are used, i.e., L/7 = 2.6 for the upper row, and L/j = 0.8 for 
the middle and lower rows. The red solid curves are computed for a true 5qp = tt/2, the blue 
dashed curves for a true 5qp = 3n/2. A true normal hierarchy is assumed. 
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energy. Comparing the beta beam and Neutrino Factory qualitatively, we find a similar behavior 
as a function of baseline and neutrino energy for both isotope pairs and all performance indicators, 
except from 18 Ne and 6 He and the 13 sensitivity. In this case, the shorter baseline dominates, which 
corresponds to the correlation only case in Fig. 3 of [38]. Quantitatively, the best performance 
for the #13 and hierarchy measurements with the Neutrino Factory comes at the magic baseline. 
The Neutrino Factory chosen here corresponds to 10 21 useful muon decays per year with muon 
energy of 50 GeV and a 50 kton magnetized iron detector. If the detector efficiency is taken as 
about 40% and background rejection factor as 1 x 10~ 5 , the sensitivity reaches in sin 2 2# 13 (true) 
for the two performance indicators #i 3 and sgr^Am 2 ^) at the magic baseline are about 2 x 10~ 4 
and 1 x 10 -4 respectively [49]. Best CP violation sensitivity comes at L ~ 3000 — 5000 km. At 
this baseline, the sin 2 2# 13 reach for maximal CP violation is 7 x 10~ 5 [49]. These numbers can be 
compared with the corresponding reach for the optimal Beta-beam set-ups identified in this paper. 
One can see that the Neutrino Factory outperforms the Beta-beam for the standard luminosity 
chosen here and for values of 7 < 650. However, for higher values of 7 and higher luminosities, 
Beta-beam returns a performance which is comparable to that of a Neutrino Factory. Recall that 
we have taken conservative estimates for the number of useful ion decays per year and the size of 
the detector. In principle, it should be possible to have at least 10 times more useful decays per 
year [61]. In addition, for some ranges of 7 it should be possible to use megaton water detectors 
with Beta-beams, and this would again increase the luminosity by factor of 10. 

6 Summary and Conclusions 

Beta-beams provide intense and well understood neutrino fluxes of a single flavor. They have 
very low beam-related backgrounds and very low systematic uncertainties, and are hence ideal 
for precision experiments. While detailed physics reach analyses of specific projects involving 
Beta-beams have been performed in the literature, greenfield set-ups have been less extensively 
studied. In this study, we have considered the Beta-beam option as a greenfield scenario and 
have identified the optimal set-ups for an iron calorimeter as the detector. We have optimized for 
maximal reaches in sin 2 2#i 3 , i.e., our Beta-beams are designed to measure sin 2 2#i 3 , CP violation, 
and the mass hierarchy for as small sin 2 2# i3 as possible. We have studied two pairs of plausible 
source ions for the Beta-beam: the standard ions 18 Ne and 6 He, which have been extensively 
studied in the literature, as well as the new candidates 8 B and 8 Li. For each pair of source 
ions, we have optimized the experimental set-up as a function of baseline L, boost factor 7, and 
luminosity. The luminosity is proportional to the number of useful ion decays, detector efficiency, 
size of detector, and exposure time. We have followed two complementary approaches for our 
study: (i) Fixing the input parameters of the beam and looking for the sensitivity reach between 
the two pairs of ions as a function of the baseline, (ii) Matching the shape (energies) of the event 
spectrum and luminosity for the two sets of ions and studying the required input parameters of 
the beam as a function of the baseline. In both approaches, we have compared the physics reach 
for the two pairs of ions in order to identify the optimal conditions for the sensitivity reach. 

In order to compare the two different pairs of isotopes, it has been useful to determine the 
conditions for the same physics output, i.e., neutrino energies, matter effects, etc., at the same 
baseline. Since the total flux at the detector for a fixed baseline is proportional to iVg7 2 (with 
N/3 being the useful ion decays per year) and the peak neutrino energy is approximately given by 
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7^0 j we have identified the conditions for obtaining similar neutrino spectral shape as N^ +Ll ~ 
12 • Np e+He , r^Ne+He ^ 35 . ^b+Li _ rj^jg means that, because of the stronger forward collimation 
of the beam, a smaller required 7 for 8 B and 8 Li has, in principle, to be compensated by a 
correspondingly larger number of useful ion decays. We have verified this behavior in realistic 
simulations. 

As the next step, we have fixed the beam input parameters Np to 1.1 x 10 18 (v e ) and 2.9 x 10 18 
(p e ) useful ion decays/year, and 7 to 350, 500, and 650, respectively. Because of the higher 
neutrino energies for the same 7, 8 B and 8 Li experience stronger matter effects than 6 He and 18 Ne. 
This has two major implications: First, if the matter resonance energy can be covered, the 1/L 2 
dependence of the events is reduced for 8 B and 8 Li, and the event rates increase at longer baselines. 
Second, because of the matter resonant or anti-resonant behavior, the mass hierarchy can be much 
easily determined. Both implications make the magic baseline attractive for the sin 2 26* 13 and mass 
hierarchy measurements using 8 B and 8 Li, has had been shown in [39, 40, 41]. For 6 He and 18 Ne, 
on the other hand, a much shorter baseline is preferable. For CP violation however, the 6 He and 
18 Ne turn out to have a discovery reach [46, 47] about one order of magnitude better than for 8 B 
and 8 Li at a relatively short baseline L ~ 600 — 1000 km (depending on the 7 chosen). Similar 
to the Neutrino Factory, we therefore observe a tension between measuring the mass hierarchy 
(preferring long baselines and high energies close to the matter resonance energy), and measuring 
CP violation (preferring short baselines determined by the oscillation maximum with as little 
matter effects as possible). Therefore, in order to optimally access all performance indicators, one 
may finally require both pairs of ions [51] and two baselines [45]. 

If one wants to compare similar neutrino energies for the two isotope pairs, i.e., one fixes the 
output, one has to include 7's about a factor of three higher for 6 He and 18 Ne than for 8 B and 
8 Li. From an optimization of our performance indicators in the L-7 plane (while still keeping the 
useful ion decays fixed), we have learned that there are, in principle, two sets of baselines which 
exhibit local optima in the performances: 

• A short baseline with L [km]/7 ~ 0.8 for 6 He and 18 Ne, and L [km]/7 ~ 2.6 for 8 B and 8 Li. 

• A long "magic" baseline L ~ 7 500 km, where the dependence on 5qp vanishes. 

While the shorter baseline is always the better choice for the CP violation measurement, the choice 
between the longer and shorter baseline for the sin 2 2# 13 sensitivity depends on 7 and the isotope 
used. For 8 B and 8 Li, the magic baseline is preferred for 7 > 350, while for 6 He and 18 Ne, the 
shorter baseline is always preferred for any realistic 7 because of the higher event rates for the 
same neutrino energies. For the mass hierarchy sensitivity, the magic baseline tends to be the best 
choice for 8 B and 8 Li already for 7 > 300, while for 6 He and 18 Ne a relatively high 7 > 1000 is 
required. However, in the latter case, the short baseline alone will not be sufficient to measure the 
mass hierarchy, another longer baseline is required. 

As the last step, we have focused on the above two sets of baselines, and we have varied the 
luminosity. For the sin 2 2#i3 sensitivity, we have demonstrated that the choice between the short 
and long baseline depends on the reference luminosity as well. For the mass hierarchy sensitivity, 
the conclusion to use a longer baseline is very robust. And for the CP violation sensitivity, which is 
only present at the short baseline, the luminosity plays a crucial role to resolve (otherwise present) 
degeneracies for 5cp = 3n/2 (for the normal hierarchy) or Sqp = vt/2 (for the inverted hierarchy). 
Note that this degeneracy resolution could also be achieved with the combination of two baselines. 
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We conclude that a greenfield Beta-beam could have excellent sensitivity reaches for the 
sin 2 2#i3, mass hierarchy, and CP violation discoveries. Comparison of the physics reach between 
( 8 B, 8 Li) and ( 6 He, 18 Ne) pairs is not at all straightforward. On the one hand, 8 B and 8 Li produce 
a given neutrino energy by a boost factor about 3.5 times lower than that needed for 18 Ne and 
6 He. Therefore, the first pair of ions would be the preferred choice if one needed higher energy 
Beta-beams within the constraint of the envisaged accelerator facilities. On the other hand, lower 
boost factors result in a lower beam collimation, and hence lower statistics, which would have to be 
compensated by a higher luminosity. For 18 Ne and 6 He one would need a 7 about 3.5 times larger 
than for 8 B and 8 Li to obtain the same neutrino energies. This constraint stretches the demanded 
7 to the prohibitively large regime. The statistics, in return, would be about a factor of 10 higher. 
Therefore, the optimal selection of ions and baselines crucially depends on the boost factor and 
luminosity used, as well as the chosen detector technology. For our reference luminosity and iron 
detector, if the Beta-beam is operated at a realistically "high" 7 ~ 500, one would optimally use 
18 Ne and 6 He at the short baseline for CP violation, 8 B and 8 Li at the magic baseline for mass 
hierarchy, and either 18 Ne and 6 He at the short baseline or 8 B and 8 Li at the magic baseline for 
sin 2 2#i3 discovery. 
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